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1   Introduction. 

The  plane  stress  solution  D-]  is  conventionally  employed 
to  estimate  the  stress  concentration  due  to  a  circular  hole  in 
a  plate  that  is  uniaxially  atretched  at  infinity.   However,  this 
solution  is  not  a  solution  of  the  exact  theory.   Nevertheless, 
we  expect  it  to  yield  an  accurate  approximation  if  e  =  ^  is 
sufficiently  "small".   We  also  expect  that  the  accuracy  will 
Increase  as  e  — j-  q. 

The  precise  relationship  between  the  plane  stress  and  the 
exact  theories  is  given  in  [j>]    for  simply  connected  plates  with 
"smooth"  boundary  curves.   In  addition,  a  boundary  layer 
procedure"'  is  outlined  for  obtaining  increasingly  accurate 
approximations  to  the  solution  of  the  exact  theory.   These 
approximations  are  given  as  "three-dimensional  corrections"  to 
the  plane  stress  solution. 

In  this  paper  we  extend  the  method  of  [3]  to  our  stress 
concentration  problem.  Results  are  obtained  in  the  form  of  a 
power  series  in  e.   We  give  here  only  terms  up  to  and  including 


We  refer  to  the  three-dimensional  linear  theory  of  elasticity 
for  homogeneous  and  isotropic  materials  as  the  exact  theory. 
There  are  some  special  cases  for  which  the  plane  stress  solution 
is  a  solution  of  the  exact  theory  [2]. 

2 

Here,  h  is  one  half  of  the  plate  thickness  and  R  is  the  radius 

of  the  hole. 

''  It  is  a  generalization  of  the  one  given  by  Friedrichs  [4]  and 
Friedrlchs  and  Dressier  [5]  in  a  study  of  the  "bending"  of  plates. 
See  also  [6]  and  [?] . 
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second  order.   Within  this  approximation  we  show  that  the  plane 
stress  theory  yields  extremely  accurate,  although  non-conservative 
predictions  of  the  maximum  stress  concentration  for  "small"  but 
finite  values  of  e.   This  accuracy  depends  upon  Polsson's  ratio  v. 
For  example,  with  v  =  I/5  the  error  in  the  plane  stress  solution 
is  less  than  5  9^)  if  e  _^  -3.   For  "larger  "  values  of  v  and  e, 
the  error  increases.   More  accurate  approximations  to  the  solution 
of  the  exact  theory,  which  may  be  necessary  for  these  values  of  e, 
can  be  obtained  by  determining  third  and  higher  terms  in  the 
expansion. 

Other  approximate  solutions  of  the  exact  theory  for  this 
problem  are  given  by  Sternberg  and  Sadowsky  [8]  using  a  modifica- 
tion of  the  Ritz  method  and  Green  [9]  and  Alblas  [10]  who 
employed  infinite  series  expansions.^  Our  results  for  the  stress 
concentration  compare  favorably  with  those  of  Alblas  if  ^  ^  jr  > 
see  Figs.  1  and  2.   Agreement  is  especially  good  with  his 
"asymptotic"  solution  which  is  closely  related  to  our  approxima- 
tion method. 

2.  Formulation. 

We  introduce  a  cylindrical  coordinate  system  r,9,z.   An 
Infinite  plate  of  thickness  2h  with  a  circular  hole  of  radius  R 
is  considered  as  a  three-dimensional  elastic  body  bounded  by  the 
planes  (the  "faces"  of  the  plate)  z  =  ±  h  and  the  cylindrical 
surface  r  =  R.   The  origin  of  the  coordinates  is  fixed  at  the 
center  of  the  hole  on  the  mldplane  of  the  plate,  z  =  0.  The  plate 
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See  also  Relssner  [11]. 
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is  stretched  at  infinity  by  a  constant  tensile  force  T  in  a 
fixed  direction  v/hich  we  take  as  the  x-axis.   The  boundary  of 
the  hole,  i.e.  the  edge,  and  the  faces  of  the  plate  are  free  of 
forces. 

If  we  introduce  the  dimensionless  variables: 


i  =  --^  ,     e  1 0  ;     c  =  f  ,     Id  1 1  , 


and  the  parameter, 

F  -  h 

then  the  faces  of  the  plate  are  given  by  i^  =  ±1  and  the  boundary 
of  the  hole  by  |  =  0. 

Considering  the  components  of  stress  as  functions  of  1,9,^ 
and  employing  an  obvious  notation,  the  stress  formulation  of  the 
exact  theory  is  given  by: 

Equilibrium  Equations, 

(1)         ^0z,c  '''^^'^re,e-^^i  +  ^^"^^^"a,9  +^^re^^=°' 


"z,C  ""^^^rzA^^^^  ^^'^^-^9z,9-'<^^rz^^    =   ^    ' 


Y-J. 


{$  ■■■■1} 


(2) 

Here, 

(3a)  A(|,e,C;e)  =-Acr  -o,^^  +  2(i+e)-^(cr^-rr9  +  2crg^9)  , 

(3b)   B(|,9,C;e)  =  -A<r0-(l+e)"^O,^+2(l+e)-2(-O,gg+2a;^-2^^-4j;g^g)  , 

(3c)  c(e,0,C;e)=  -AT:0-(i+l)"^o,^g+(i+e)"^(o,g+4cr3+a^^g-f);^g)  , 
Acr=  a-,^.+(i+e)-^ot|+(i+U"^^,0e  >   ^  =  TT^^^r"^  ^e +^)  ^ 

V  Is  Poisson's  ratio  and  a  comma  indicates  partial  differentiatior 
To  complete  the  formulation  we  require  appropriate  boundary 
conditions.   These  are  obtained  by  specifying  the  applied  forces 

as :  „ ,  ,  ,.  .   ,  , 

(4a)   cr^^(l,9,±l)  =  a^^(^,9,±l)=  .r^(!;,0,±l)  =  0  ; 

(4b)   (T  ^(0,9,0=  :^(O,0,a  =  org(O,9,a  =  0  ; 

(4c)  <r^^{^,e,a   =0    ,      .X,(-,9,C)=Tcos^9,3^g(co,9,0  =-  ^  T  sin  29  , 

where  T  is  the  constant  tensile  force  at  infinity  in  the  x- 
direction. 


.  t 
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Since  we  are  concerned  only  with  the  extension  of  plates, 
it  can  be  shown  [5],  v/ithout  loss  of  generality,  that  cf  ,    t'q,  cr 
and  CTg  are  even  functions  of  C  while  ^      and  ^  are  odd 
functions  of  C'   In  the  following  sections  we  shall  make  frequent 
use  of  these  properties  without  explicit  reference. 

3.  The  Interior  Problems. 

We  assume  that  each  stress  component,  indicated  by  the 
generic  symbol  a'(|,9,^;£),  can  be  asymptotically  represented  by 
a  power  series  in  e: 

CO 

(5)         -^(1,0, C;b)  -I~r"(?,9,C)£''  . 

n^ 

Here  <r     are  called  the  interior  stress  coefficients  and  we  define 
a^  =  0  if  n  <  0.   We  assume  that  the  jr    possess  the  same  even  or 
odd  property  as  cr.   Substituting  these  expansions  into  (1-5) 
and  (4a)  and  equating  coefficients  of  like  powers  of  e  yields  a 
system  of  differential  equations  that  are  satisfied  by  the  ^y^ . 
The  analysis  of  this  system  is  elementary  and  similar  to  that 
outlined  in  [3]-   Therefore,  the  calculations  are  not  explicitly 
exhibited.   Instead  some  of  the  results  are  listed  below.   For 
example,  we  can  show  that: 


(7) 


-l^^   +  (i+e)-^(.rje,e-^-^?  - 'A^  =  "^  >  V  ^=0,1, 


(ft-i^ 


Is.f.c*n-3T»' 


•  Ctfg&'j 


(8a)   (r"  =  s"(?,0),  cr"}  =  S^{^,6),     cr-"  =  S^.(^,0),   n=0,l; 


r9 


^n  ^  1  ^n-2(^^gj^2^3n(^^g)^  ^^n  ^  1  3n-2(  ^^^^  ^2  ^  gH^  ^^^  j 


(8b) 


where  a"^,  b"^  and  C^  are  obtained  from  (3)  using  (5). 

We  define  tlie  m-th  order  Interior  problem  (Problem  I  )  as 
the  boundary  value  problem  that  contains  (6-S)  with  n  =  m  j^  3 
and  appropriate  boundary  conditions  which  are  obtained  In  the 
following  sections.   Equations  (6-8)  with  n  =  0  are,  in  our 
notation,  the  stress  relations  and  the  differential  equations 
of  the  classical  theory  of  plane  stress  in  polar  coordinates  [1]. 

k.    Formulation  of  the  Boundary  Layer  Problem. 

The  results  given  in  the  previous  section  are  obtained 
without  reference  to  the  edge  boundary  conditions  (4b).   In  fact, 
the  expansions  (5)  cannot,  in  general,  satisfy  these  conditions. 
If  they  did,  it  then  follows  from  (8)  and  (3)  that  more  boundary 
conditions  then  are  appropriate  for  the  solution  of  (7)  are 
specified.   Thus,  if  the  series  (5)  represent  the  solution  of 
the  exact  theory  they  do  so  in  a  region  away  from  the  edge  which 
we  call  the  "interior  domain".   The  region  of  the  plate  adjacent 
to  and  including  the  edge  where  (5)  may  deviate  rapidly  from  the 
solution  of  the  exact  theory  is  called  the  "boundary  layer".   To 
obtain  expansions  that  may  converge  uniformly  up  to  and  including 
the  edge  we  assume,  as  in  [3]  ,  that  the  deviation  of  (5)  from 

See  also  [6,7] • 
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the  exact  solution  occurs  only  In  the  direction  normal  to  the 
edge,  i.e.  in  the  ^  direction.   V/e  then  introduce  the  "stretched" 
boundary  layer  variable  [12]  t]   as: 

(9)  r,=l. 

Considering  the  r],d,t^   coordinate  system,  boundary  layer 
stresses  indicated  by  the  generic  symbol  f(T),9,^;e)  are  defined 
as: 

(10)  f(ri,9,(:;e)  =<T'{e,9,C;e)  • 

To  determine  approximations  to  the  exact  stress  distribution 
near  the  edge  of  the  hole  and  appropriate  boundary  conditions 
for  the  interior  problems  we  assume  that: 


(11)  f(Ti,9,C;£)  ^YZ   f''(n,9,03"  . 

n-0 

The  f  are  called  the  boundary  layer  stress  coefficients  and  we 

n  n 

set  f  =  0  if  n  *=  0.   We  further  assume  that  each  f"  possesses 

the  same  even  or  odd  property  as  the  corresponding  f. 

Substituting  (9-10)  into  the  exact  theory  (l-'+)  and  equating 

coefficients  of  like  powers  of  e  we  find  from  the  coefficients 

of  e"  that  the  f"  satisfy: 

f?z,/:  ^f?,n  +IZ'(fJe,e  -^  fj;  -  f^)  =  o  , 


viudi-'- 


.r.irc-   "'-f"'",'-'  noi^tyloE   .-trpr"^  srfJ 


'b-rir. 


W-).      ~  V 


'.■;.■  7,r:.L       "■ 


^rul^L   ^"1.3 


sbht/on 


v2^n.  nn 


II  r    . 


r    '    'lit)     *■ —     V ,T\     ^—    L  r,90  r      9        r9,9 


=   0    , 


(13) 


,71  'Tj 


r      0        r0,9'    '•9,90    I     '09_ 


=   0    , 


^2  n  ^pn  J-'f^         +>~"(f^      ^a-f^   -2fi      ^)    =   0    , 

rz    '     'ri^       ^-—     rz,ri     ' —        rz,  99      rz        9z,9'  ' 


72^n 


-^"9  +  11    (fr9,,  +  nl,^9)-i:    Ls^9,99+2^(^r--f9)-^fr9-nl9j    =   0, 


Jz,Ti"i     '9C' 


,4,99-1z-^24,9l    =    0    > 


(14) 


f^^(0,9,a    =   f^(0,9,C)    =   f2g(O,0,C)    =   0    , 


Cjoo,9,a    =   0    , 


fS(co,9,a    = 


T  cos   9,      n  =  0 
0  ,      n  ^  0 


Here, 


r 


„  (••  -T/2   sin  29,      n  =   0 

^0  ,      n  ?^  0 


_1_  (fH  n  nv  ^2  _     ^^  5^ 


1       T t /    ,  vJ^JaI  t—  a^ 


i+j+l=n 


YZ-^Tl  (-i)^{j+i)tiJ'a=^  . 

i+j-i?=n 


In  addition  we  require  the  boundary  layer  stresses  to 
approach  or  "match"  the  interior  stresses  as  e  — >  o. 
Specifically,  we  assume  that  each  o''^(4>9,0  has,  near  ^  =  0, 
a  Taylor  series  expansion  in  |: 


I^  t. 


to 


-  \^ 


2t 


i?i;^lJfr 


where 


(15a)  3n,  j.^B-'^lO  S,0  n=0,l,...  . 

Therefore  from  (5)  and  (9)  we  have,  In  some  region  about  ^  =  0, 

00  ^ 

(16)  cr{^,0.^;e)  ^rir''(T,9,Oe''  , 

n=0 

where 

(15b)  ^''(^,9,0  =  V—  s;;-"^(9,an'" 

m=0 

are  the  Interior  coefficients  near  ^  =  0  as  functions  r\ ,    9   and  C* 
For  each  n  we  define  the  ''redu'-'ed  boundary  layer  stress 
coefficients",  F^(r],9,0  as: 

(17)  P'^(n,e,a  s  f''(n,e,0  --^""(n^e^C)  . 

The  "matching  condition"  or  the  asymptotic  form  for  the  f  is 
obtained  from  (9),  (10)  and  (l6)  by  associating  each  f"  with  the 

corresponding  o^  as  e  — >  0.   Using  (9)  and  ( ?.? )  we  write  this 

i 
condition  as: ' 

(15c)  lim  F^(t],9,C)  =  0,  n=0,l,.... 

Tj — >C0 

An  analysis,  simi?-c^r  to  the  preceding,  applied  to  the 
boundary  at  "infinity"  yields  a  corresponding  formulation  of  that 
boundary  layer  problem.   We  do  not  exhibit  this  formulation. 

'  In  obtaining  (15g)  we  assume  that  the  terms  in  f"  which  \anish 
a-s  T]   — *.  00  do  so  faster  than  any  negative  power  of  r) . 


( ■'.  v.^)'^-, 


%:1±d: 
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5.  Analysis  of  the  Boundary  Layer  Problems. 

For  each  n  Eqs.  (12-15)  and  (1?)  separate  into  two 

distinct  systems  which  we  call  Problem  P^  and  Problem  T  . 

Problem  P'^  is  concerned  with  the  coefficients  f^,  f"^,  f2  and  f" 

r   z   y      rz 

and  involves  the  first  two  of  (12),  the  first  four  of  (13),  (15) 
and  (1?)  for  these  coefficients  and  the  first  two  of  each  of  (14). 
Problem  T  is  concerned  with  the  remaining  two  coefficients  f  g 
and  Tq     and  the  remaining  equations  in  (12-15)  and  (17). 

We  shall  associate  with  p'^  a  "stress  function",  ^   (il,9,C), 
which  may  be  the  solution  of  the  following  boundary  value  problem 
on  the  sem.i-infinite  strip,  |^1  _^  1,  t^  >_  0,  and  fixed  9: 

V^f  =  0  ; 


i?^^(n,^±l)  =  <t^?.,^^(il,0,±l)  =  0;     iim   [(!)?.,  ^?^^]  =   0  ; 

<!)':, .(o,e,c)  =  P'^(e,c)  ,  (t>';'„^{o,e,c)  =  o  , 


where  6"'S^<^)  =  p  (0,-C).  It  can  be  shovm  by  ^.lerr.entary  means 
.,  .  -  .m    .n. 


thrt  if  6';r.^,  I',;  p  and  h''^\,   are  uniformly  con'.iouous  functions 
of  C  and  if  I  ,  ^,      and  i,^  are  single-valued  functions  then. 


(19)  f    p''^(9,(;)ciC  =  0  for  all  9    . 


Eiriployirg  Problei.i  l\    we  can  show  that  a  solution  of  P  is 
given  by 

(20a;       F_  =  9, pp.    F   =  6,   ,    F    =  -Q ,    f, ,      Fa  =  vv  (h 


OH'S  3     (TI 

S'H  m'^L<io-i'i  bap 


i  0  -^  -f^ 


:        1 .'  J 


=  iii.'y^r 


;  bar   o      ^    . 
bi-.s    ./?(|;    ^''} 


pr ) 
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If  (j)  is  the  solution  of  (l8)  with  n  =  0  and 
(20b)  P°(e,a  =  -S°(0,9)  . 

A.  0 

The  first  boundary  condition  on  ^  =  0  for  Problem  I   is  obtained 


from  (20b)  and  (19)  as. 


(21)  S°(0,9)  =  0  . 


Using  Problem  I   it  follows  that 


is  a  solution  of  Problem  T  If  f   {t]  ,9 ,1^)   ±s   the  solution  of. 


(23) 


V  ^  =  0   , 


^?^(Ti,e,±i)  =  0,   lim  M.,   ^^  ]  =  0,   ^?  (o,9,c)  =g"(e,0 


with  n  =  0  where. 


(22b)  g°(e,C)  =  -s°0(o,0)  . 


Single  valuedness  of  the  solution  of  (23)  yields  from  (22b)  the 
second  boundary  condition  on  |  =  0  for  I  as, 

(24)  S°g(O,0)  =  0  . 

From  (18)  and  (20-24)  it  follows  that 

(25)  F°(Ti,e,C)  =  0  . 


0  -   %^7 


.0  - 


0  -   ( 
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In  a  similar  manner  it  may  be  shown  from  an  analysis  of 
Problems  P"^  and  T  that 

(26)  F^(Ti,e,0  5  0 

and  that  the  boundary  conditions  on  ?  =  0  for  Problem  I   are 
given  by, 

(27)  Sl{0,9)    =   sIq{0,9)    =  0  . 

A  solution  of  Problem  P^  is  obtained  from  (l8)  with  n  =  2 
and 
(28a)  P^(9,a  =  -S2(0,e).i  aO(0,9)c2 

if 

(28b)   F^-^?,,,  ^1-^%,'     ^l.  =  <X'     ^e  =  ^^'^'- 

Equations  (19)  and  (28a)  yield  the  first  boundary  condition  for 

2 
Problem  I  as, 

(29)  S2(0,9)  =  -  I  A°(O,0) 

where  A°  is  obtained  from  (5)  and  (:5a).   Similarly  we  obtain  a 
solution  of  Problem  T^  by  setting 

O  P        2      2 

Pre  =  -  ^'c  '  ^ez  =  ^'n 

where  ^^(ti,9,C)  is  a  solution  of  (25)  with  n  =  2  and 


'O    CiriO  £•'.  X 


G'j 


■  r     ,'  "'      '^        '^ 
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Single  valuedness  of  the  solution  gives  the  second  boundary 

2 
condition  for  Problem  I  as. 


(30) 


S;q{0,6)   = 


-  -i  c°(o,e) 


More  accurate  approximations  of  the  stresses  near  the  edge  can  be 
obtained  by  examining  Problems  P  and  T  for  n  ^  J>, 

A  corresponding  analysis  of  the  boundary  layer  at  "infinity' 
yields  the  appropriate  boundary  conditions  at  infinity  for  the 
interior  problems.   For  example,  we  can  show  that, 

(31a)    S°(oo,e)  =  |(1+  cos  29)  ,    S^g( ro , 0)  =  -  |  sin  20  , 


(31b)    S^(co,e)  =  S^g(co,0)  =  0 


if   n  =  1,2  . 


6 .  Solution  of  the  Interior  and  Bounda ry  Layer  Problems . 

Problem  I   (the  plane  stress  theory)  which  consists  of  the 
differential  equations  (7)  and  (8a)  with  n  =  0  and  the  boundary 
conditions  (21),  (24)  and  (31a)  has  the  solution  [l]: 


2  0  _  2  „0  _  T 
T  ^r  -  T  ^r  "  ^  - 


+ri+-^ 


(32)    ^Tq   =  -^   Sg  =  l  +  __^.  [_. 


(1+1)^   ^    (1+^)^^   (l+D^J 
^— -jf  cos  29 


cos  20  , 


1  + 


2  _0   _  2  „0 
T"  r0  "  T  ^r0 


^   (1+ 


(1+4) 
—2^  sin  20  , 


e)"   (1+4) 


rz    9z    z 


Grid    'it^"!    \, 


yS  file 
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Since  the  differential  equations  (7)  and  the  boundary  conditions 
(27)  and  (31b)  are  homogeneous  we  have  for  the  solution  of 


Problem  I  : 


(33)  0-1  =  0. 


2 
Employing  (32)  and  (3)  we  obtain  the  solution  of  Problem  I  , 

which  is  given  by  (7)  and  (Sb)  with  n  =  2  and  (29),  (30)  and 
(31b),  as 


(34a) 


where 


G-^  =  -  GTg  =  G  cos  29,   a-^9  =  G  sin  29, 


2     2     2   ^ 


(34b)  G{^,u  -jYm-jj:;:^  ' 

2 
To  obtain  a  solution  of  Problem  P  we  must  solve  the 

2 
boundary  value  problem  (l8)  with  n  =  2  where  ^  is  given  by 

(2Ca),  (29),  (3a)  and  (32).   Since  an  exact  solution  of  this 

problem  is  unknown  we  employ  the  "approximate"  solution  given 

by  Horvay  [13] : 

.  ,■  .  .  ■       ■  ,.  •,       ,  '  ■  •' 

p2  =  ^e'^'^lcosbri  +-|  sin  bri)(3C^-l)  cos  29  , 

2   2  •  •  <■ 

f2  =  H(^^-)  e"^"^  (-cos  bn+^  sin  bTi)(l-C^)2  cos  29  , 
z      ic;  D 

(35a) 

2   2 

f2  =  H(5-i^)  e"^"^  sin  bo  (C  -C^)  cos  29  , 
rz       3d 


F9  =  ^(^?  +  ^z^  > 


where 


fcrfs   I  Of)    <(Q5^)    r  :ip  ■   brus   (T)   vcf  fi9vx3  si  dolrivgF 


,es  soo  0  =  |x  -  =  ^- 


n    ^   =>  ^^  -■' 


(Br^) 


9"e3riv; 


v^-iT  "  (3.^)0  {ci.^O 


$:ij   •  a^Xdoil  lo  no.WwXoa  b  ni:B>>cfo  oT 

vd  n^vj"  .-•-..-fr-  !>  -  ri  rfctxw  {8l}  mslcfcriq  sjjXbv  •'^^fim/ocf 

iioviig  CiOJ:.-  •  .•  smxxo'iqqB "  srid-  YoXqnrs  sw  iiwoinlny  el  nieXdCTg 

,   bS  aoo   (l-^:S^)(r.d  nl8  |+ prd  soo)'^®"3  p  =  ^"5 

5  /  S  -      ,    ,  ,  S   ,        ,  s     fTB-        /     d->-    iJ ,  tr  Srr 
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(55b)       H  =  1^  ,   a  =  2.075  ,  b  =  1.143  . 


An  Integral  representation  can  be  given  for  the  solution 

2 
of  Problem  T  since  Green's  function  for  the  boundary  value 

problem  (23)  is  knovm  [l4].   However,  we  prefer  to  use  the 

Infinite  series  representation  which  yields  for  the  solution 
.2  ' 


of  T' 


(36a) 


where. 


rt  oo  p 


P^^  =  >_  K     cos  nTT?   e"'"'^   sin  29      , 
re       ^     n 


^Q     =  YH  K^   sin  nTrC   e"'^^^   sin  29      , 


(36b)  ,-      •         K     =  ^«   (-1)" 


n        ^2  2 

TT         n 


As  in  [3]  we  define  the  stresses,  r       ,  of  the  N-th 
approximation  to  the  exact  theorj'-  as: 


(37)    (T^^h^,9,^;E)   =rZ[s^{^^,9,r)  +  F^I,^/e,e,0]e''    • 


n=0 

.(0) 


It  follows  from  this  definition,  (25),  (26)  and  (33)  that  j 
and  (T         coincide  with  the  plane  stress  solution.   For  the 
second  approximation  to  the  exact  theory  we  obtain  from  (37) 


fioX^if J.O&  •s-dv'-  -^col  ns-zis  sd  hbo  I'ioiii  .  -^  Ib':J2 

Q 

•  ,'oa  ^fi!f  Tol  aftXf'l'r  rio-tdv?  norti?:':  v:  ssi'^foe  9.-J-J-r:.!:'>rf-?. 


,     feS  iix3  ^^""s  ?-trn  eoo  ^^^I  2» 


V'  .1 


[b'K] 


^^'V     ■:t£rf,1    i^^).  fcn  .(eS)    ^rioictirrnsfo  sln^  mc'.;r  c:;.'Oj..L 01  jI 

(L) 


(38) 
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cr^^^4,e,C;e)  =  s°(4, 9)  +  [.71(4,9,0  +  F^(^/e,9,0]e^  , 
4i^4,e,^;e)  =  sO0(4,9)  +  [cT^e(4,9,C)  +  P^0(?/e,9,O]e^  , 
3'i2)(4,9,C;e)  =  F^  {e/£,9,ae^  , 


4z^(^>^-J£)  =  F^^(|/e,e,C)e^  , 


0-1^^(4, 9, C;e)  =  f2(4/£,0,O£^  . 


,0    2    ,  „2 


Here  S  ,  Gr    and  F  are  given  in  (32)  and  (3^-36). 


7.  Presentation  of  Results. 

We  define  the  quantity: 


(39)      ^^^^4,C;£)  =  ^ 5^^^^ 


1  +-^ 


(1+?J'^. 


In  Fig.  1,  Oq    '(0,0;e)  is  illustrated  as  a  function  of  e  with 

varying  Poisson's  ratio.   For  e  =  0  we  obtain  the  plane  stress 

^(2) 

result  which  gives  a  stress-concentration  factor  of  two  for  cfg 

Although  the  stress-concentration  factor  increases  with  e,  for 
"small"  e,  it  is  only  a  small  percentage  of  the  plane  stress 
solution.   For  example,  with  e  =  -2  and  v  =  1/4  our  results 
Indicate  only  a  1.25  9^  increase  over  the  plane  stress  result. 
Thus,  for  "small"  e  the  plane  stress  solution  apparently  yields 


i^r 


,  ^sriJ.G-^-^^'^- 


.-  1 ..  .: 


(2) 


'T> 


,  ''3[(3\-e<3\'^)e^'5+(?.e,^)^--cj 


ml 


.  rtvVi.w-.- 


y/fg-  s-iB  ^"^  fonj5  S:>  t'*8  s>'59H 


j'lj^  j- V 


--^ 


fi  8B  be:f£iieu£Li 


.•blues'!  eeisncte  sfislq  srict  'lavo  aafisioni    c/?  .Hii! 


no  61   if  I    <3      ilBrna'' 
")  to"!!     ;noi:tt/Xo8 
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accurate  but  unconservative  results.   The  circles  in  Pig-  1 
represent  the  results  of  Alblas  [10]   (v  =  1/4)  obtained  from 
a  formal  infinite  series  solution  of  the  exact  theory.   The 
dotted  curve  gives  the  "asymptotic  solution"  of  Alblas. 

Figure  2  reveals  the  behavior  of  cr^   (0,±1;£)  as  a 
function  of  e.   The  curve  obtained  from  the  "asymptotic  solution" 

of  Alblas  (v  =  1/4  and  e  ^   1/4)  coincides  with  our  curve.   In 

^  (o) 

Fig.  3  the  variation  through  the  thickness  of  ^^    '    is  indicated 

at  the  edge  of  the  hole  for  v  =  1/4  and  varying  e.   Since  cTg 
assumes  its  maximum  on  the  middle  surface,  these  results  are  of 
some  importance  in  experiments  where  measurements  are  taken  on 
the  faces  of  the  plate. 

In  the  remaining  three  figures  we  illustrate  "boundary 
layer  behaviors"  for  some  of  the  stresses.   In  Fig.  4  the 
variation  with  e,  of  the  middle  surface  values  of  cTq        is  given 
for  V  =  1/4  and  two  values  of  s.   The  l-variation  of  3"^  '  is 

relatively  insensitive  to  changes  in  e  especially  for  e  >    .2. 

'■■ '  0 ) 
Using  the  plane  stress  solution  we  define  T'g   as 


^0    1 

(40)  j(o).ZJ^ 


1  + 


(1+^)^^ 


cos"  20 

This  quantity  is  in  close  agreement  with  Tq    hi,0;e)    for  e  =  l/lO 
and  is  not  shown  in  Fig.  4.        ,,,... 

Figure  5  shows  the  ^-variation  of  the  scaled  thickness 
shear  stress, 


b 
1C    dTTS    eilU^OT    9831".^-      ^  :K<JiSS 


^■■•T  fyhcWcd 


Oi\.^ 


:woria    jura  Ki  dhb 
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(41)         S-^^Ut  C's)  -( 1.99jLll+yXr ')  rf^^hi  9  C'e) 

^^^'  ^z9    ^^'^'^f    -\    2kvT   sin  29  )  ^z9    ^^'^'^»^' 

for  C  =  1/2.   We  observe  that  for  e  =  l/k,     ^^^y    at  C  =  .1  is 
only  J>0  9o  of  its  maximum  value,  while  at  ?  =  '2,  (T  g'  is 
8.59^  of  its  maximum  value.   If  we  define  the  "boundary  layer 
thickness  ,  ^  ,  as  that  value  of  s  at  which  CT  q     is  a  snail 
percentage,  say  5  9^  >  of  its  maximum  value  then  the  results 
illustrated  in  Pig.  5  give 

In  the  remaining  graph.  Fig.  6,  the  ^   variation  of  <^  a     is 


.(2) 
r9 


defined  as, 

.(2) 


ry  /_     T  sm  26    ^   (1+1)^   (1+4)^1  j 


/v  f  2 ) 

For  the  plane  stress  solution,  e  =  0,  the  corresponding  o'  a 

coincides  with  the  |  axis. 

We  wish  to  emphasize  that  more  accurate  approximations ,to 
the  exact  solution  can  be  obtained  by  extending  our  calculations 
to  terms  of  order  three  and  greater.   Additional  accuracy  in  our 

solution  (38)  could  be  obtained  if  the  exact  solution  to  Problem 

2 
P  ,  rather  than  (35),  were  available.   Since  (36)  provides  the 

2  (2)        (2) 

exact  solution  to  Problem  T  ,  the  stresses  ^\q     and  CTg  '  are 

most  likely  given  with  greater  precision  than  the  other  stresses 
in  (38). 


3?   *■ 


(^o.®<-      '  ')4lTi         txn 


lo   n'"'.ri'^ 


^    |.^r-^^   -  ■  -{^J)^ 
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2.20  r 


2.15 


See  Eq.  (39) 


2.00 


Figure  1:   Variation  uith  e  and  v  of  the  raiddle  plane 
values  of  &t;' '    at  the  edge  of  the  hole. 


22 


2.0 


1.9 


See  Eq.  (39) 


V  -  lA 


Figure  2:   Variaticn  v;ith  e  and  v  of  the  face  plane 
values  of  ^]y      at  the  edge  of  the  hole. 
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See  Eq.  (39) 


V  =  lA 


s  =  1/10 
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Fi.  ure  3:   Variation  through  the  thickness  of 


.(2) 
9 

V   =   lA. 


CTA""    at    the    edge    of    the    hole    for 


C'!^ 
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See  Eqs.  (39)  and  (1+0) 


V  =  lA 


Figure  l\.:      Variation  i-ilth   i   of  the  middle  plane  values 
(2) 


of  ^Q^    for  V  =  lA. 
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See    jilq.    (ij.1) 


Figure   S-      Variation  with   ^ 


of  a^Q^  foi"  ^  =  1/2 
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See  Eq.  (i|2) 


V  =  lA 


!      \         \    e  =  1/3 


e  =  lA 


"^\  e  =  lAo 


Figure  6:   Variation  with  £,   of  a^^^  for  ^  =  1, 
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Washington  7,   D.  C. 

Attn:   Code  l40 

(1) 

Code  600 

(1) 

Code  700 

(1) 

Code  720 

(1) 

Code  725 

(1) 

Code  731 

(1) 

Code  740 

(2) 

CO,  U.S.  Naval  Ordnance  Lab. 
White  Oak,  Maryland 
Attn:   Techn.  Library  (2) 

Techn.  Eval.  Dept.       (1) 

Director,  Materials  Lab. 

N.Y.  Naval  Shipyard 

Brooklyn  1,  N.  Y.  (1) 

CO,  Port.smouth  Naval  Shipyard 
Portsmouth,  New  Hampshire       (2) 

CO,  Mare  Island  Nav.  Shipyard 
Vallejo,  California  (2) 

CO  and  Director 

U.S.  Nav.  Electron.  Lab. 

San  Diego  52,  California        (1) 

Officer-ln-Charge 
Nav,  Civ.  Engin.  Res. 
e.nd  Eval.  Lab. 
U.S.  Nav.  Constr.  Battal.  Center 
Port  Hu.rneme,  California        (2) 

Dir. ,  Nav.  Air  Experimental  Sta. 
Nav.  Air  Mat.  Center,  Nav.  Base 
Philadelphia  12,  Penna. 
Attn:   Materials  Lab.  (1) 

Structures  Lab.  (l) 

Officer-ln-Charge 

Underwater  Explos.  Res.  Dlv. 

Norfolk  Naval  Shipyard 

Portsmouth,  Virginia 

Attn:   Dr.  A.  H.  Kell  (2) 

CO,  U.S.  Nav.  Proving  Ground 
Dahlgren,  Virginia  (1) 

Supr.  of  Shipbuilding 

USN  and  Nav.  Inspec.  of  Ordnance 

General  Dynamics  Corp. , 

Electr.  Boat  Dlv. 

Grcton,  Connecticut  (1) 

Supr.  of  Shipbuilding 

USN  and  Nav.  Inspec.  of  Ordnance 

Newport  News  Shipbuilding 

and  Dry  Dock  Co. 
Newport  News,  Virginia         (1) 
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Supr.  of  Shipbuilding 
USN  and  Nav.  Inspec.  of  Ordnance 
Ingalls  Shipbuilding  Corp. 
Pascagoula,  Mississippi      (1) 

CO,  U.S.  Nav.  Admin.  Unit 

MIT,  Cambridge  59,  Mass.     (1) 

Officer-ln-Charge 
Postgrad.  School  for 

Naval  Officers 
Webb  Inst,  of  Nav.  Arch. 
Crescent  Beach  Rd. 
Glen  Cove,  L.I.,  N.Y.        (l) 

Supt.,  Nav.  Gun  Factory 
Washington  25,  D,  C.         (1) 

Comm. ,  Nav.  Ordnance  Test  Sta. 
China  Lake,  California 
Attn:   Physics  Div.  (l) 

Mechanics  Div.        (1) 

CO,  Nav.  Ordnance  Test  Sta. 

Underwater  Ordnance  Div. 

3202  E.  Foothill  Blvd. 

Pasadena  8,  California 

Attn:   Struc.  Div.  (1) 

CO  and  Director 

U.S.  Nav.  Engin.  Exp.  Station 

Annapolis,  Maryland 

Supt.  U.S.  Nav.  Postgrad 
Monterey,  California 

Comm.  Marine  Corps  Schools 

Quantico,  Virginia 

Attn:   Dlr. ,  MC  Dev.  Center   (l) 

Comm.  Gen. ,  USAF 
Washington  25,  D.  C. 
Attn:   Res.  and  Dev 


(1) 

School 
(1) 


Div. 

CO,  Air  Material  Command 
Wright-Patterson  APB,  Ohio 
Attn:   MCREX-B 

Structures  Div. 

CO,  USAF  Inst,  of  Technology 
Wright-Patterson  AFB,  Ohio 
Attn:   Chief,  Appl.  Mech. 
Group 


(1) 


(1) 
(1) 


(1) 


Director  of  Intelligence 

Headquarters,  USAF 

Washington  25,  D.  C. 

Attn:   PV  Br.  (Air  Targ.  Div)(l) 

CO,  A.F.  Office  Sci.  Research 

Washington  25,  D.  C. 

Attn:   Mechanics  Div.        (1) 


U.S.  Atomic  Energy  Commission 
VJashington  25,  D.  C. 
Attn:   Dir.  of  Research 

Dir.,  Nat.  Bur.  of  Standards 
Washington  25,  D.  C. 
Attn:   Div.  of  Mechanics 

Engin.  Mech.  Sect. 

Aircraft  Structures 

Comm.,  U.S.  Coast  Guard 

1300  E  St. ,  NW 

Washington  25,  D.  C. 

Attn:   Chief,  Test  and  Dev.  Div 

U.S.  Maritime  Administration 
General  Admin.  Office  Bldg. 
441  G  St.,  NW 
Washington  25,  D.  C. 
Attn:   Cnief,  Div.  Prelim. 
Design 

Nat.  Aero,  and  Space  Admin. 
Langley  Research  Center 
Langley  Field,  Virginia 
Attn:   Structures  Div. 

Nat.  Aero,  and  Space  Admin. 

1512  H  St, ,  NW 

Washington  25,  D.  C. 

Attn:   Loads  and  Struc.  Div. 

Director,  Forest  Prod.  Lab. 
Madison,  VJisconsin 

Federal  Aviation  Agency 
Dept.  of  Commerce 
Washington  25,  D.  C. 
Attn:   Chief,  Air  Engin.  Div. 

Chief,  Air.  and 

Equip.  Div. 

National  Science  Foundation 
1520  H  St. ,  m 
Washington,  D.  C, 

National  Academy  of  Sciences 

2101  Constitution  Ave., 

Washington  25,  D.  C. 

Attn:   Dir.,  Comm.  on  Ships 
Struc.  Design 
Exec.  Secy.,  Comm.  on 
Undersea  Warfare 

General  Dynamics  Corp. 
Electr.  Boat  Div. 
Groton,  Connecticut 

Newport  News  Shipbuilding  and 

Dry  Dock  Co. 
Newport  News,  Virginia 
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(1) 
(1) 
(1) 


(1) 


(1) 
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(2) 
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(1) 
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(1) 

(1) 
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Ingalls  Shipbuilding  Corp. 
Pascagoula,  Mississippi      (1) 

Prof.  Lj-nn  S.  Beedle 

Fritz  Engineering  Lab. 

Lehigh  University 

Bethlehem,  Penna.  (1) 

Prof.  R.  L.  Bisplinghoff 
Dept.  of  Aero.  Engineering 
Massachusetts  Inst,  of  Techn. 
Cambridge  39,  Massachusetts   (1) 

Prof.  H.  H.  Blelch 
Dept.  of  Civ.  Engineering 
Columbia  University- 
New  York  27,  N.  Y.  (1) 

Prof.  B.  A.  Boley 

Dept.  of  Civ.  Engineering 

Columbia  University 

New  York  2?,  N.  Y.  (l) 

Dr.  John  P.  Brahtz 

Southern  California  Labs. 

Stanford  Research  Institute 

820  Mission  St. 

South  Pasadena,  California   (l) 

Dr.  D.  0.  Brush 

Struc.  Dept.  53-13 

Lockheed  Aircraft  Corp. 

Missile  Syst.  Div. 

Sunnyvale,  California        (l) 

Prof.  B.  Budiansky 

Dept.  of  Mech.  Engineering 

School  Appl.  Sciences 

Harvard  University 

Cambridge  38,  Massachusetts   (1) 

Prof.  Herbert  Deresiewicz 

Dept.  of  Civ.  Engineering 

Columbia  University 

632  W.  125th  St. 

New  York  27,  N.  Y.  (1) 

Prof.  D.  C.  Drucker,  Chmn. 

Div.  of  Engineering 

Brown  University 

Providence  12,  Rhode  Island   (1) 

Prof.  John  Duberg 

Dept.  of  Civ.  Engineering 

University  of  Illinois 

Urbana,  Illinois  (l) 

Prof.  J.  Ericksen 
Mech.  Engineering  Dept. 
Johns  Hopkins  University 
Baltimore  I8,  Maryland       (1) 


Prof.  A.  C.  Eringen 

Dept.  Aero.  Engineering 

Purdue  University 

Lafayette,  Indiana  (l) 

Prof.  W.  Flugge 

Dept.  of  Mech.  Engineering 

Stanford,  California  (1) 

Mr.  M.  Goland,  VP  and  Dir. 
Southwest  Research  Institute 
8500  Culebra  Rd. 
San  Antonio  6,  Texas  (l) 

Prof.  J.  N.  Goodier 

Dept.  of  Mech.  Engineering 

Stanford  University 

Stanford,  California  (1) 

Prof.  L.  E.  Goodman 
Engineering  Experimental  Sta. 
University  of  Minnesota 
Minneapolis,   Minnesota         (1) 

Prof.  M.  Hetenyi 

The  Technical  Institute 

Northwestern  University 

Evanston,  Illinois  (1) 

Prof.  P.  G.  Hodge 

Dept.  of  Mechanics 

Illinois  Inst,  of  Technology 

Chicago  16,  Illinois  (1) 

Prof.  N,  J.  Hoff,  Head 

Div.  Aeronautical  Engineering 

Stanford  University 

Stanford,  California  (1) 

Prof.  W.  H.  Hoppmann,  II 

Dept.  of  Mechanics 

Rensselaer  Polytechnic  Inst. 

Troy,  New  York  (l) 

Prof.  Bruce  G.  Johnston 

University  of  Michigan 

Ann  Arbor,  Michigan  (l) 

Prof.  J.  Kempner 
Dept.  of  Aero,  Engineering 
and  Appl.  Mechanics 
Polytechnic  Inst,  of  Brooklyn 
333  Jay  St. 
Brooklyn  1,  N.  Y.  (1) 

Prof.  H.  L.  Langhaar 
Dept.  of  Theoretical  and 
Applied  Mechanics 
University  of  Illinois 
Urbana,  Illinois  (1) 
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Prof.  B.  J.  Lazan,  Director 
Engineering  Experimental  Sts. 
University  of  Minnesota 
Minneapolis  l4,  Minnesota    (1) 

Prof.  E.  H.  Lee 

Dlv.  of  Appl.  Mathematics 

Brovm  University 

Providence  12,  Rhode  Island   (1) 

Prof.  George  H.  Lee,  Dlr.  of  Res. 
Rensselaer  Polytechnic  Inst. 
Troy,  New  York  (l) 

Mr.  S.  Levy 

GE  Electr.  Research  Lab. 

3198  Chestnut  St. 

Philadelphia,  Penna.         (1) 

Prof.  Paul  Lleber 
Geology  Department 
University  of  California 
Berkeley  4,  California       (l) 

Prof.  Joseph  Marin,  Head 
Dept.  Engineering  Mechanics 
College  of  Engln.  and  Arch. 
Pennsylvania  State  University 
University  Pai^k,  Penna.      (l) 

Prof.  R.  D.  Mindlin 

Dept.  of  Civ.  Engineering 

Columbia  University 

632  W.  125th  St. 

New  York  2?,  N.  Y.  (l) 

Prof.  Paul  M.  Naghdl 

Building  T~7 

College  of  Engineering 

University  of  California 

Berkeley  4,  California       (l) 

Prof.  William  A.  Nash 
Dept.  of  Engineering  Mechanics 
University  of  Florida 
Gainesville,  Florida         (l) 

Prof.  N.  M.  Newmark,  Head 
Dept.  of  Civ.  Engineering 
University  of  Illinois 
Urbana,  Illinois  (l) 

Prof.  E.  Orowan 
Dept.  of  Mech.  Engineering 
Massachusetts  Institute  of  Techn. 
Cambridge  39,  Massachusetts   (1) 

Prof.  Aris  Phillips 

Dept.  of  Civ.  Engineering 

15  Prospect  St. 

Yale  University 

New  Haven,  Connecticut       (l) 


Prof.  W.  Prager,  Chmn. 

Phys.  Sol.  Council 

Brovm  University 

Providence  12,  Rhode  Island     (1) 

Prof.  J.  R.  M.  Radok 
Dept.  of  Aero  Engineering 

and  Appl.  Mechanics 
Polytechnic  Inst,  of  Brooklyn 
333  Jay  St. 
Brooklyn  1,  N.  Y.  (1) 

Prof.  E.  L.  Relss 

Inst,  of  Mathematical  Sciences 

New  York  University 

4  Washington  Place 

New  York  3,  N.  Y.  (1) 

Prof.  E.  Reissner 
Dept.  of  Mathematics 
Massachusetts  Inst,  of  Technology 
Cambridge  39,  Massachusetts     (l) 

Prof.  M.  A.  Sadowsky 

Dept.  of  Mechanics 

Rennselaer  Polytechnic  Inst. 

Troy,  New  York  ( l) 

Prof.  B.  W.  Shaffer 

Dept.  of  Mech.  Engineering 

New  York  University 

University  Heights 

New  York  53,  N.  Y.  (1) 

Prof.  J.  Stallmeyer 

Dept.  of  Civ.  Engineering 

University  of  Illinois 

Urbana,  Illinois  (1) 

Prof.  Eli  Sternberg 

Dept.  of  Mechanics 

Brown  University 

Providence  12,  Rhode  Island     (1) 

Prof.  T.  Y.  Thomas 

Grad.  Inst.  Math,  and  Mech. 

Indiana  University 

Bloomlngton,  Indiana  (1) 

Prof.  S.  P.  Timoshenko 

School  of  Engineering 

Stanford  University 

Stanford,  California  (1) 

Prof.  A.  S.  Velestos 

Dept.  of  Civ.  Engineering 

University  of  Illinois 

Urbana,  Illinois  (1) 
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Dr.  E.  Wenk 

Southwest  Research  Institute 

8500  Culebra  Rd. 

San  Antonio,  Texas  (1) 

Prof.  Dana  Young 
Yale  University- 
New  Haven,  Connecticut       (l) 

Prof.  R.  A.  Di  Taranto 

Dept.  of  Mech.  Engineering 

Drexel  Institute 

52nd  and  Chestnut  Streets 

Philadelphia,  Penna.         (l) 

Mr.  H.  K.  Koopman,  Secy. 
Welding  Res.  Council 

Engineering  FounCation 
29  W.  39th  St. 
New  York  I8,  N.  Y.  (2) 

Prof.  Walter  T.  Daniels 
School  of  Engin.  and  Archit. 
Howard  University 
VJashington  1,  D.  C.  (1) 

Comm. ,  (Code  755) 

U.S.  Naval  Ordnance  Test  Sta. 

China  Lake,  California 

Attn:   Techn.  Library        (1) 

Prof.  J.  E.  Cermak 

Dept.  of  Civ.  Engineering 

Colorado  State  University 

Fort  Collins,  Colorado       (1) 

Prof.  W.  J.  Hall 

Dept,  of  Civ.  Engineering 

University  of  Illinois 

Urbana,  Illinois  (1) 

Dr.  Hyman  Serbin 

Design  Integration  Dept. 

Hughes  Aircraft  Co. 

Culver  City,  California      (1) 

Commander 

Wright  Air  Development  Center 

Wright-Patterson  Air  Force  Base 

Dayton,  Ohio 

Attn:  Dynamics  Branch  (l) 
Aircraft  Lab.  (1) 
WCLSY  ( 1 ) 

Commanding  Officer 

USNNOEU 

Kirtland  Air  Force  Base 

Albuquerque,  New  Mexico 

Attn:   Code  20 

(Dr.  J.N.  Brennan)    (l) 


Legislative  Reference  Service 
Library  of  Congress 
Washington  25,  D.  C, 
Attn:   Dr.  E.  Wenk 


1) 


Dr.  A.  Ross 

Aircraft  Nuclear  Propulsion  Dept. 

General  Electric  Co. 

Cincinnati  I5,  Ohio  (1) 

Dr.  P.  Lane 

General  Applied  Science  Labs. 
Stewart  and  Merrick  Avenues 
Westbury,  L.I.,  N.  Y.  (1) 

Commander  Edv/ard  Leonard 

Asst.  Navy  Representative 

MIT  Lincoln  Lab. 

Lexington  75>  Massachusetts     (1) 
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